      SUBROUTINE  UVXYAI(C,C2,XI,M,JEN,JIN,NN,LL,UE1,UE3,UI1,VE1,VE3,   00010000
     +                   VI1,XE1,XE3,XI1,YE1,YE3,YI1,IPO)               00020006
CCC           INTEGRALS A-I  &  PARAMETERS U,V,X,Y                      00030000
CCC           PARAMETERS U,V,X &Y  -  COMPLEX                           00040000
CCC           IPO = 1   PROLATE SYSTEM                                  00050000
CCC           IPO =-1   OBLATE SYSTEM                                   00060000
CCC    ----------ORIGINALLY PROGRAMMED IN 1973 ; REVISED IN 1977------- 00070000                                          00070000
      PARAMETER  (JMX=150,NMX= 70)                                      00071010
      COMPLEX*16 UE1,UE3,UI1,VE1,VE3,VI1,XE1,XE3,XI1,YE1,YE3,YI1,AI1,   00080006
     1           AI2,AI3,AI4,AI5,AI6,AI7,AI8,AI9,DIN1,DIN2,DIN3,DIN4,   00090006
     2           DIN5,DIP1,DIP3,DIP5,DI(NMX,JMX),RE3(NMX),DRE3(NMX),    00100009
     3           RI1(NMX),DRI1(NMX),SPCI(NMX),C1,C2,CI,CJ               00110005
      REAL*8     AE1,AE2,AE3,AE4,AE5,AE6,AE7,AE8,AE9,DEN1,DEN2,DEN3,    00120008
     1           DEN4,DEN5,DEP1,DEP3,DEP5,XI,XI2,DE(NMX,JMX),           00130009
     2           SPCE(NMX),RE1(NMX),DRE1(NMX),C,FLOAT                   00140008
      COMMON  /XYUV/SPCI,RI1,DRI1,RE3,DRE3,DE,RE1,DRE1,SPCE             00150005
      COMMON  /CMPLX/DI                                                 00151005
      DATA   CI,CJ/(0.D0,1.D0),(1.D0,0.D0)/                             00160000
C                                                                       00170000
      FLOAT(I)=DFLOAT(I)                                                00180000
      C1=CJ*C                                                           00190000
      XI2=XI*XI-1.D0                                                    00200000
       IF(IPO.EQ.-1) XI2=XI*XI+1.D0                                     00210003
      NM=(M+NN-1)-M                                                     00220000
      I=LL-1                                                            00230000
      IM1=I+M+M                                                         00240000
      IM2=I+IM1                                                         00250000
      II=I                                                              00260000
       IF(MOD(I,2).EQ.0) II=I+1                                         00270003
C             INTEGRAL  A - I                                           00280000
      AE1 = 0.D0                                                        00290000
      AE2 = 0.D0                                                        00300000
      AE3 = 0.D0                                                        00310000
      AE4 = 0.D0                                                        00320000
      AE5 = 0.D0                                                        00330000
      AE6 = 0.D0                                                        00340000
      AE7 = 0.D0                                                        00350000
      AE8 = 0.D0                                                        00360000
      AE9 = 0.D0                                                        00370000
      AI1 = (0.D0,0.D0)                                                 00380000
      AI2 = (0.D0,0.D0)                                                 00390000
      AI3 = (0.D0,0.D0)                                                 00400000
      AI4 = (0.D0,0.D0)                                                 00410000
      AI5 = (0.D0,0.D0)                                                 00420000
      AI6 = (0.D0,0.D0)                                                 00430000
      AI7 = (0.D0,0.D0)                                                 00440000
      AI8 = (0.D0,0.D0)                                                 00450000
      AI9 = (0.D0,0.D0)                                                 00460000
          IF(MOD(NM+I,2).EQ.0)  THEN                                    00470003
        IF(M.GT.0)  THEN                                                00480003
      DEN2=0.                                                           00490000
      DEN4=0.                                                           00500000
      DIN2=(0. ,0. )                                                    00510000
      DIN4=(0. ,0. )                                                    00520000
       IF(I-2.LT.0) GO TO 15                                            00530003
      DEN2=DE(NN,II-2)                                                  00540000
      DIN2=DI(NN,II-2)                                                  00550000
       IF(I-4.LT.0) GO TO 15                                            00560003
      DEN4=DE(NN,II-4)                                                  00570000
      DIN4=DI(NN,II-4)                                                  00580000
   15 AE1=FLOAT((IM1-1)*IM1)*DE(NN,II)/FLOAT(IM2+1)-FLOAT(I*(I-1))*     00590000
     +    DEN2/FLOAT(IM2-3)                                             00600005
      AI1=FLOAT((IM1-1)*IM1)*DI(NN,II)/FLOAT(IM2+1)-FLOAT(I*(I-1))*     00610000
     +    DIN2/FLOAT(IM2-3)                                             00620005
      AE2=FLOAT(IM2-1)*DE(NN,II)                                        00630000
      AI2=FLOAT(IM2-1)*DI(NN,II)                                        00640000
      AE6=-FLOAT(I*(I-1)*(IM2+1)+IM1*(IM1-1))*DE(NN,II)/FLOAT(2*(IM2+1))00650000
     +    -FLOAT(I*(I-1)*(I+M-2))*DEN2/FLOAT(IM2-3)                     00660005
      AI6=-FLOAT(I*(I-1)*(IM2+1)+IM1*(IM1-1))*DI(NN,II)/FLOAT(2*(IM2+1))00670000
     +    -FLOAT(I*(I-1)*(I+M-2))*DIN2/FLOAT(IM2-3)                     00680005
      III=II+2                                                          00690000
      DO 2000 KA=III,JEN,2                                              00700000
      AE2=AE2+FLOAT(IM2-1)*DE(NN,KA)                                    00710000
 2000 AE6=AE6+FLOAT(M*(IM2-1))*DE(NN,KA)                                00720000
      DO 2010 KB=III,JIN,2                                              00730000
      AI2=AI2+FLOAT(IM2-1)*DI(NN,KB)                                    00740000
 2010 AI6=AI6+FLOAT(M*(IM2-1))*DI(NN,KB)                                00750000
      AE7=FLOAT((IM1-1)* IM1*(IM1+1)*(IM1+2))/FLOAT((IM2+1)*(IM2+3))*   00760000
     +   (DE(NN,II)/FLOAT(IM2+1)-DE(NN,II+2)/FLOAT(IM2+5))-2.0*FLOAT(I* 00770005
     +   (I-1)*(IM1-1)*IM1)/FLOAT((IM2-3)*(IM2+1))*(DEN2/FLOAT(IM2-3)-  00780005
     +   DE(NN,II)/FLOAT(IM2+1))+FLOAT((I-3)*(I-2)*(I-1)*I)/FLOAT(      00790005
     +   (IM2-3)*(IM2-5))*(DEN4/FLOAT(IM2-7)-DEN2/FLOAT(IM2-3))         00800005
      AI7=FLOAT((IM1-1)* IM1*(IM1+1)*(IM1+2))/FLOAT((IM2+1)*(IM2+3))*   00810000
     +   (DI(NN,II)/FLOAT(IM2+1)-DI(NN,II+2)/FLOAT(IM2+5))-2.0*FLOAT(I* 00820005
     +   (I-1)*(IM1-1)*IM1)/FLOAT((IM2-3)*(IM2+1))*(DIN2/FLOAT(IM2-3)-  00830005
     +   DI(NN,II)/FLOAT(IM2+1))+FLOAT((I-3)*(I-2)*(I-1)*I)/FLOAT(      00840005
     +   (IM2-3)*(IM2-5))*(DIN4/FLOAT(IM2-7)-DIN2/FLOAT(IM2-3))         00850005
        END IF                                                          00860003
          ELSE                                                          00870003
        IF(MOD(I,2).EQ.0)  THEN                                         00880003
      IK=I+1                                                            00890000
      DEP1=DE(NN,I+1)                                                   00900000
      DIP1=DI(NN,I+1)                                                   00910000
      DEP3=DE(NN,I+3)                                                   00920000
      DIP3=DI(NN,I+3)                                                   00930000
      DEP5=DE(NN,I+5)                                                   00940000
      DIP5=DI(NN,I+5)                                                   00950000
      DEN1=0.                                                           00960000
      DEN3=0.                                                           00970000
      DEN5=0.                                                           00980000
      DIN1=(0. ,0. )                                                    00990000
      DIN3=(0. ,0. )                                                    01000000
      DIN5=(0. ,0. )                                                    01010000
       IF(I-1.LT.0) GO TO 10                                            01020003
      DEN1=DE(NN,I-1)                                                   01030000
      DIN1=DI(NN,I-1)                                                   01040000
       IF(I-3.LT.0) GO TO 10                                            01050003
      DEN3=DE(NN,I-3)                                                   01060000
      DIN3=DI(NN,I-3)                                                   01070000
       IF(I-5.LT.0) GO TO 10                                            01080003
      DEN5=DE(NN,I-5)                                                   01090000
      DIN5=DI(NN,I-5)                                                   01100000
   10  CONTINUE                                                         01110003
        ELSE                                                            01120003
      IK=I+2                                                            01130003
      DEP1=DE(NN,I+2)                                                   01140000
      DIP1=DI(NN,I+2)                                                   01150000
      DEP3=DE(NN,I+4)                                                   01160000
      DIP3=DI(NN,I+4)                                                   01170000
      DEP5=DE(NN,I+6)                                                   01180000
      DIP5=DI(NN,I+6)                                                   01190000
      DEN1=0.                                                           01200000
      DEN3=0.                                                           01210000
      DEN5=0.                                                           01220000
      DIN1=(0. ,0. )                                                    01230000
      DIN3=(0. ,0. )                                                    01240000
      DIN5=(0. ,0. )                                                    01250000
       IF(I-1.LT.0) GO TO 20                                            01260003
      DEN1=DE(NN,I)                                                     01270000
      DIN1=DI(NN,I)                                                     01280000
       IF(I-3.LT.0) GO TO 20                                            01290003
      DEN3=DE(NN,I-2)                                                   01300000
      DIN3=DI(NN,I-2)                                                   01310000
       IF(I-5.LT.0) GO TO 20                                            01320003
      DEN5=DE(NN,I-4)                                                   01330000
      DIN5=DI(NN,I-4)                                                   01340000
   20  CONTINUE                                                         01350003
        END IF                                                          01360003
         IF(M.GT.0)  THEN                                               01370003
      AE3=FLOAT(IM1*(IM1-1))*(FLOAT(IM1+1)*DEP1/FLOAT(IM2+3)+FLOAT(I)*  01380000
     +    DEN1/FLOAT(IM2-1))/FLOAT(IM2+1)-FLOAT(I*(I-1))*(FLOAT(IM1-1)* 01390005
     +    DEN1/FLOAT(IM2-1)+FLOAT(I-2)*DEN3/FLOAT(IM2-5))/FLOAT(IM2-3)  01400005
      AI3=FLOAT(IM1*(IM1-1))*(FLOAT(IM1+1)*DIP1/FLOAT(IM2+3)+FLOAT(I)*  01410000
     +    DIN1/FLOAT(IM2-1))/FLOAT(IM2+1)-FLOAT(I*(I-1))*(FLOAT(IM1-1)* 01420005
     +    DIN1/FLOAT(IM2-1)+FLOAT(I-2)*DIN3/FLOAT(IM2-5))/FLOAT(IM2-3)  01430005
      AE4=FLOAT(I)*DEN1                                                 01440000
      AI4=FLOAT(I)*DIN1                                                 01450000
      AE5=-FLOAT(I*(I+M-1))*DEN1                                        01460000
      AI5=-FLOAT(I*(I+M-1))*DIN1                                        01470000
      DO 2020 KC=IK,JEN,2                                               01480000
      AE4=AE4+FLOAT(IM2-1)*DE(NN,KC)                                    01490000
 2020 AE5=AE5+FLOAT(M*(IM2-1))*DE(NN,KC)                                01500000
      DO 2030 KD=IK,JIN,2                                               01510000
      AI4=AI4+FLOAT(IM2-1)*DI(NN,KD)                                    01520000
 2030 AI5=AI5+FLOAT(M*(IM2-1))*DI(NN,KD)                                01530000
      AE8=FLOAT((IM1-1)*IM1*(IM1+1)*(IM1+2)*(IM1+3))/FLOAT((IM2+1)*     01540000
     +  (IM2+3)*(IM2+5))*(DEP1/FLOAT(IM2+3)-DEP3/FLOAT(IM2+7))+FLOAT(I* 01550005
     +  (M+M-I)*(IM1-1)*IM1*(IM1+1))/FLOAT((IM2-3)*(IM2+1)*(IM2+3))*(   01560005
     +  DEN1/FLOAT(IM2-1)-DEP1/FLOAT(IM2+3))-FLOAT((I+4*M-1)*(IM1-1)*I* 01570005
     +  (I-1)*(I-2))/FLOAT((IM2-5)*(IM2-3)*(IM2+1))*(DEN3/FLOAT(IM2-5)- 01580005
     +  DEN1/FLOAT(IM2-1))+FLOAT(I*(I-1)*(I-2)*(I-3)*(I-4))/FLOAT((IM2- 01590005
     +  3)*(IM2-5)*(IM2-7))*(DEN5/FLOAT(IM2-9)-DEN3/FLOAT(IM2-5))       01600005
      AI8=FLOAT((IM1-1)*IM1*(IM1+1)*(IM1+2)*(IM1+3))/FLOAT((IM2+1)*     01610000
     *  (IM2+3)*(IM2+5))*(DIP1/FLOAT(IM2+3)-DIP3/FLOAT(IM2+7))+FLOAT(I* 01620005
     *  (M+M-I)*(IM1-1)*IM1*(IM1+1))/FLOAT((IM2-3)*(IM2+1)*(IM2+3))*(   01630005
     *  DIN1/FLOAT(IM2-1)-DIP1/FLOAT(IM2+3))-FLOAT((I+4*M-1)*(IM1-1)*I* 01640005
     *  (I-1)*(I-2))/FLOAT((IM2-5)*(IM2-3)*(IM2+1))*(DIN3/FLOAT(IM2-5)- 01650005
     *  DIN1/FLOAT(IM2-1))+FLOAT(I*(I-1)*(I-2)*(I-3)*(I-4))/FLOAT((IM2- 01660005
     *  3)*(IM2-5)*(IM2-7))*(DIN5/FLOAT(IM2-9)-DIN3/FLOAT(IM2-5))       01670005
      AE9=FLOAT(I*(I-1))/FLOAT(IM2-3)*(FLOAT((I-2)*(I+M-3))*DEN3/FLOAT  01680000
     +   (IM2-5)-FLOAT((I+M)*(IM1-1))*DEN1/FLOAT(IM2-1))-FLOAT((IM1-1)* 01690005
     +   IM1 )/FLOAT(IM2+1)*(FLOAT(I*(I+M-1))*DEN1/FLOAT(IM2-1)-FLOAT(  01700005
     +   (I+M+2)*(IM1+1))*DEP1/FLOAT(IM2+3))                            01710005
      AI9=FLOAT(I*(I-1))/FLOAT(IM2-3)*(FLOAT((I-2)*(I+M-3))*DIN3/FLOAT  01720000
     *   (IM2-5)-FLOAT((I+M)*(IM1-1))*DIN1/FLOAT(IM2-1))-FLOAT((IM1-1)* 01730005
     *   IM1 )/FLOAT(IM2+1)*(FLOAT(I*(I+M-1))*DIN1/FLOAT(IM2-1)-FLOAT(  01740005
     *   (I+M+2)*(IM1+1))*DIP1/FLOAT(IM2+3))                            01750005
         ELSE                                                           01760003
      AE3=(FLOAT(I+1)*DEP1/FLOAT(2*I+3)+FLOAT(I)*DEN1/FLOAT(I+I-1))/    01770003
     +    FLOAT(I+I+1)-(FLOAT(I+3)*DEP3/FLOAT(I+I+7)+FLOAT(I+2)*DEP1/   01780006
     +    FLOAT(I+I+3))/FLOAT(I+I+5)                                    01790006
      AI3=(FLOAT(I+1)*DIP1/FLOAT(2*I+3)+FLOAT(I)*DIN1/FLOAT(I+I-1))/    01800000
     +    FLOAT(I+I+1)-(FLOAT(I+3)*DIP3/FLOAT(I+I+7)+FLOAT(I+2)*DIP1/   01810006
     +    FLOAT(I+I+3))/FLOAT(I+I+5)                                    01820006
      AE5=DEP1                                                          01830000
      AI5=DIP1                                                          01840000
      K=I+I                                                             01850000
      AE8=FLOAT((I+3)*(I+4)*(I+5))/FLOAT((K+5)*(K+7))*(DEP1/FLOAT((K+3) 01860000
     +  *(K+5))-2.0*DEP3/FLOAT((K+5)*(K+9))+DEP5/FLOAT((K+9)*(K+11)))+  01870007
     +  3.0*FLOAT(I*(I+3))/FLOAT((K+1)*(K+5))*(DEN1/FLOAT((K-1)*(K+1))  01880007
     +  -2.0*DEP1/FLOAT((K+1)*(K+5))+DEP3/FLOAT((K+5)*(K+7)))-FLOAT(I*  01890007
     +  (I-1)*(I-2))/FLOAT((K-1)*(K+1))*(DEN3/FLOAT((K-5)*(K-3))-2.0*   01900007
     +   DEN1/FLOAT((K-3)*(K+1))+DEP1/FLOAT((K+1)*(K+3)))               01910007
      AI8=FLOAT((I+3)*(I+4)*(I+5))/FLOAT((K+5)*(K+7))*(DIP1/FLOAT((K+3) 01920000
     +  *(K+5))-2.0*DIP3/FLOAT((K+5)*(K+9))+DIP5/FLOAT((K+9)*(K+11)))+  01930007
     +  3.0*FLOAT(I*(I+3))/FLOAT((K+1)*(K+5))*(DIN1/FLOAT((K-1)*(K+1))  01940007
     +  -2.0*DIP1/FLOAT((K+1)*(K+5))+DIP3/FLOAT((K+5)*(K+7)))-FLOAT(I*  01950007
     +  (I-1)*(I-2))/FLOAT((K-1)*(K+1))*(DIN3/FLOAT((K-5)*(K-3))-2.0*   01960007
     +   DIN1/FLOAT((K-3)*(K+1))+DIP1/FLOAT((K+1)*(K+3)))               01970007
      AE9=(FLOAT((I+1)*(I+2))*DEP1/FLOAT(K+3)-FLOAT(I*(I-1))*DEN1/FLOAT 01980000
     *    (K-1))/FLOAT(K+1)-(FLOAT((I+3)*(I+4))*DEP3/FLOAT(K+7)-FLOAT(  01990007
     *    (I+1)*(I+2))*DEP1/FLOAT(K+3))/FLOAT(K+5)                      02000007
      AI9=(FLOAT((I+1)*(I+2))*DIP1/FLOAT(K+3)-FLOAT(I*(I-1))*DIN1/FLOAT 02010000
     *    (K-1))/FLOAT(K+1)-(FLOAT((I+3)*(I+4))*DIP3/FLOAT(K+7)-FLOAT(  02020007
     *    (I+1)*(I+2))*DIP1/FLOAT(K+3))/FLOAT(K+5)                      02030007
         END IF                                                         02040003
          END IF                                                        02050003
C             PARAMETERS U,V,X AND Y                                    02060000
      UE1 = (0.D0,0.D0)                                                 02070000
      UE3 = (0.D0,0.D0)                                                 02080000
      UI1 = (0.D0,0.D0)                                                 02090000
      VE1 = (0.D0,0.D0)                                                 02100000
      VE3 = (0.D0,0.D0)                                                 02110000
      VI1 = (0.D0,0.D0)                                                 02120000
      XE1 = (0.D0,0.D0)                                                 02130000
      XE3 = (0.D0,0.D0)                                                 02140000
      XI1 = (0.D0,0.D0)                                                 02150000
      YE1 = (0.D0,0.D0)                                                 02160000
      YE3 = (0.D0,0.D0)                                                 02170000
      YI1 = (0.D0,0.D0)                                                 02180000
        IF(IPO.EQ. 1) THEN                                              02190003
         IF(MOD(NM+I,2).EQ.0) THEN                                      02200003
        IF(M.GT.0) THEN                                                 02210003
      UE1=CJ*FLOAT(M)*XI*RE1(NN)*(XI2*XI2*AE2+2.0*XI2*AE1+AE7)          02220000
      UE3=CJ*FLOAT(M)*XI*RE3(NN)*(XI2*XI2*AE2+2.0*XI2*AE1+AE7)          02230000
      UI1=CJ*FLOAT(M)*XI*RI1(NN)*(XI2*XI2*AI2+2.0*XI2*AI1+AI7)          02240000
      YE1=CI*FLOAT(M)*(RE1(NN)*(AE1+AE6)/XI2+(RE1(NN)+XI*DRE1(NN))*AE2) 02250000
     +    /C1                                                           02260007
      YE3=CI*FLOAT(M)*(RE3(NN)*(AE1+AE6)/XI2+(RE3(NN)+XI*DRE3(NN))*AE2) 02270000
     +   /C1                                                            02280007
      YI1=CI*FLOAT(M)*(RI1(NN)*(AI1+AI6)/XI2+(RI1(NN)+XI*DRI1(NN))*AI2) 02290000
     +   /C2                                                            02300007
        END IF                                                          02310003
         ELSE                                                           02320003
      VE1=CI*(FLOAT(M*M)*RE1(NN)*(XI2*XI2*AE4+2.0*XI2*AE3+AE8)/XI2-     02330003
     *    SPCE(NN)*RE1(NN)*(XI2*AE3+AE8)+XI*XI2*DRE1(NN)*(2.0*AE3+XI2*  02340007
     *    AE5+AE9)+RE1(NN)*(XI2*XI2*AE5+(3.0*XI*XI-1.0)*AE9))/C1        02350007
      VE3=CI*(FLOAT(M*M)*RE3(NN)*(XI2*XI2*AE4+2.0*XI2*AE3+AE8)/XI2-     02360000
     *    SPCE(NN)*RE3(NN)*(XI2*AE3+AE8)+XI*XI2*DRE3(NN)*(2.0*AE3+XI2*  02370007
     *    AE5+AE9)+RE3(NN)*(XI2*XI2*AE5+(3.0*XI*XI-1.0)*AE9))/C1        02380007
      VI1=CI*(FLOAT(M*M)*RI1(NN)*(XI2*XI2*AI4+2.0*XI2*AI3+AI8)/XI2-     02390000
     *    SPCI(NN)*RI1(NN)*(XI2*AI3+AI8)+XI*XI2*DRI1(NN)*(2.0*AI3+XI2*  02400007
     *    AI5+AI9)+RI1(NN)*(XI2*XI2*AI5+(3.0*XI*XI-1.0)*AI9))/C2        02410007
      XE1=CJ*(XI*RE1(NN)*AE5-DRE1(NN)*AE3)                              02420000
      XE3=CJ*(XI*RE3(NN)*AE5-DRE3(NN)*AE3)                              02430000
      XI1=CJ*(XI*RI1(NN)*AI5-DRI1(NN)*AI3)                              02440000
         END IF                                                         02450003
        ELSE                                                            02460003
         IF(MOD(NM+I,2).EQ.0)  THEN                                     02470003
        IF(M.GT.0)  THEN                                                02480003
      UE1=CJ*FLOAT(M)*XI*RE1(NN)*(XI2*XI2*AE2-2.0*XI2*AE1+AE7)          02490000
      UE3=CJ*FLOAT(M)*XI*RE3(NN)*(XI2*XI2*AE2-2.0*XI2*AE1+AE7)          02500000
      UI1=CJ*FLOAT(M)*XI*RI1(NN)*(XI2*XI2*AI2-2.0*XI2*AI1+AI7)          02510000
      YE1=CI*FLOAT(M)*(AE2*(RE1(NN)+XI*DRE1(NN))-RE1(NN)*(AE1+AE6)/XI2) 02520000
     +    /C1                                                           02530007
      YE3=CI*FLOAT(M)*(AE2*(RE3(NN)+XI*DRE3(NN))-RE3(NN)*(AE1+AE6)/XI2) 02540000
     +    /C1                                                           02550007
      YI1=CI*FLOAT(M)*(AI2*(RI1(NN)+XI*DRI1(NN))-RI1(NN)*(AI1+AI6)/XI2) 02560000
     +    /C2                                                           02570007
        END IF                                                          02580003
       ELSE                                                             02590003
      VE1=CI*(SPCE(NN)*RE1(NN)*(XI2*AE3-AE8)-FLOAT(M*M)*RE1(NN)*        02600003
     *   (XI2*XI2*AE4-2.0*XI2*AE3+AE8)/XI2+XI*XI2*DRE1(NN)*(XI2*AE5-AE9 02610007
     *   -2.0*AE3)+RE1(NN)*(XI2*XI2*AE5-(3.0*XI*XI+1.0)*AE9))/C1        02620007
      VE3=CI*(SPCE(NN)*RE3(NN)*(XI2*AE3-AE8)-FLOAT(M*M)*RE3(NN)*        02630000
     *   (XI2*XI2*AE4-2.0*XI2*AE3+AE8)/XI2+XI*XI2*DRE3(NN)*(XI2*AE5-AE9 02640007
     *   -2.0*AE3)+RE3(NN)*(XI2*XI2*AE5-(3.0*XI*XI+1.0)*AE9))/C1        02650007
      VI1=CI*(SPCI(NN)*RI1(NN)*(XI2*AI3-AI8)-FLOAT(M*M)*RI1(NN)*        02660000
     *   (XI2*XI2*AI4-2.0*XI2*AI3+AI8)/XI2+XI*XI2*DRI1(NN)*(XI2*AI5-AI9 02670007
     *   -2.0*AI3)+RI1(NN)*(XI2*XI2*AI5-(3.0*XI*XI+1.0)*AI9))/C2        02680007
      XE1=CJ*(XI*RE1(NN)*AE5+DRE1(NN)*AE3)                              02690000
      XE3=CJ*(XI*RE3(NN)*AE5+DRE3(NN)*AE3)                              02700000
      XI1=CJ*(XI*RI1(NN)*AI5+DRI1(NN)*AI3)                              02710000
         END IF                                                         02720003
        END IF                                                          02730003
      RETURN                                                            02740000
      END                                                               02750000
